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We review the theoretical ideas and tools required to arrive at a next-to-next-to- 
leading order (NNLO) description of the heavy quark-antiquark production cross 
section in e+e — annihilation for the case that the center of mass kinetic energy 
of the quarks is larger than Aqcd- In this case the NNLO cross section can be 
calculated with purely perturbative methods. We present details of the calculation 
and discuss two applications, the determination of the bottom quark mass from T 
sum rules and the top-antitop production cross section close to threshold in e+e — 
annihilation. 

1 Introduction 

Within the last year significant progress has been achieved in the perturbative 
description of heavy quark-antiquark (QQ) pairs in the kinematic regime close 
to threshold using the concepts of effective field theories. In this talk we report 
on some of these developments from the point of view of practical applications. 
There have also been many interesting new results concerning the systematics 
of the effective field theory description of QQ pairs at threshold which would 
deserve to be presented in detail but can only be mentioned peripherally due 
to lack of space. 

To be definite, we consider the total cross section a thr = a(e + e~ — » QQ + 
anything) for the cm. energies yfs ss 2Mq, Mq being the heavy quark mass. 
In this kinematic regime the QQ dynamics is characterized by the fact that the 
heavy quarks have small cm. velocities. Considering the strong dependence 



of the velocity (v — y 1 — AMq/s) in the threshold regime on Mq it is obvious 
that, at least in principle, accurate calculations of the cross section may lead to 
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precise determinations of the quark mass. The situation, however, is not simple 
due to the fact that the QQ system in the nonrelativistic regime is governed by 
(at least) three scales. This makes this system particularly complicated. The 
three scales are Mq, the quark mass, Mqv, the relative momentum, and Mqv 2 , 
the kinetic energy in the cm. frame. This has three important consequences 
which are symptomatic for all nonrelativistic QQ systems: 

1. Because the three scales are widely separated (Mq ^> Mqv » 
Mqv 2 ). The theoretical tools required to describe the QQ pair strongly de- 
pend on whether the scales each are smaller or larger then Aqcd- The most 
transparent situation arises if Mqv 2 > Aqcd because in this case the theoreti- 
cal tools resemble most closely those of QED bound state calculations. In fact, 
for the two applications mentioned in this talk it is ensured that this condition 
is satisfied. 

2. Because ratios of the three scales arise, the description of the nonrel- 
ativistic QQ pair involves a double expansion in a s and v. This means that 
the standard multi-loop expansion in a s breaks down. The most prominent 
indication of this fact is the so called Coulomb singularity, originating from 
the ratio Mq/Mqv, which corresponds to a singular (a s /t;)™-behavior in the 
n-loop correction to the amplitude 7 — > QQ for v — > 0. The latter singular- 
ity is caused by the 00-component of the gluon propagator and, in Coulomb 
gauge, is directly related to the exchange of longitudinally polarized gluons. 
This singularity has to be treated by resumming diagrams involving longitudi- 
nal gluon exchange to all orders in a s . The Coulomb singularity exists in the 
nonrelativistic limit, but there are also power- like and logarithmic divergences 
in v which are suppressed by powers of the velocity. At this point we would like 
to specify more clearly what "resummation" means in this context. Strictly 
speaking it would mean a resummation of the perturbation series in a s to all 
orders, where the respective coefficients are expanded up to a certain power 
in v. This means that the resummation would be carried out in the (formal) 
limit a s « d « 1. The resulting series would then (uniquely) define analytic 
functions which could be continued to the region of interest \v\ <~ a s . Typical 
structures like QQ bound states can only be observed after this continuation. 
Of course, this would be a highly cumbersome and inefficient method. It is 
therefore mandatory to reformulate the problem in terms of wave equations. 
The solutions of these wave equations are equivalent to the results of the re- 
summation method. As a matter of convenience we will also call the wave 
equation method "resummation" for the rest of this talk. 

3. From point 1 we can conclude that the cm. velocity of the heavy quarks 
should satisfy the condition v > (Aqcd /Mq) 1 / 2 to ensure that the scale Mqv 2 
is perturbative. In addition, relativistic corrections are not suppressed by fac- 
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tors of 7r like multi-loop corrections. This means that relativistic corrections 
can be quite sizeable. The corrections of 0(v 2 ), called NNLO from now on, 
can be estimated to be of order 20 — 30% for bb and 5% for ti. Thus, the calcu- 
lation of higher order relativistic corrections is mandatory in order to achieve 
sufficient theoretical accuracy and to test the reliability of the perturbative 
description itself. Obviously the perturbative treatment works better if Mq is 
large. 

In order to arrive at a reliable theoretical description of a thr we have to go 
through two steps: first, we have to address the question how to organize the 
calculation systematically keeping in mind points 1-3, and, second, we have to 
actually carry out the calculation itself. In Section ^ we will briefly address 
the systematics and in Section || we will present the calculation of the photon 
mediated cross section at NNLO in the nonrelativistic expansion. Section [| is 
devoted to two applications, the determination of the bottom quark mass and 
the ti cross section at threshold in e + e~ annihilation. 



2 Systematics and NRQCD 

A very economical approach to systematically deal with the problems described 
previously is to take advantage of the separation of the scales Mq, Mqv and 
Mqv 2 using the concepts of effective field theories. In the following we outline 
the conceptual steps to arrive at a NNLO description of a nonrelativistic QQ 
pair for Aqcd < Mqv 2 without going very far into formal considerations. The 
basic idea of the effective field theory approach is to integrate out momenta 
above the scales relevant for the nonrelativistic dynamics of the QQ pair. Doing 
this, one always has to keep in mind the relation of each of the scales to Aqcd ■ 
(If Mq were of order as or even smaller than Aqcd a nonrelativistic expansion 
would be meaningless in the first place.) 

Suppose that Mq is larger than Aqcd- In this case we can integrate out 
momenta of order Mq because they are not responsible for the nonrelativistic 
dynamics of the QQ pair. Starting from QCD we then arrive at an effective 
field theory in which the heavy quarks and the gluons interacting with them 
only carry momenta below Mq . This forces us to introduce different fields for 
heavy quark and antiquark. The resulting theory is called nonrelativstic QCD 
(NRQCD) El and its Lagrangian reads 
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The gluonic and light quark degrees of freedom are described by the con- 
ventional relativistic Lagrangian, whereas the heavy quark and antiquark are 
described by the Pauli spinors tp and x> respectively. For convenience all color 
indices are suppressed. Only those terms relevant for the NNLO cross section 
are displayed, where we have omitted the straightforward antitop bilinears. 
The latter can be obtained through charge conjugation symmetry. The effects 
coming from momenta or order Mq are encoded in the short-distance coeffi- 
cients Oi, . . . ,05. They can be determined as a perturbative series in a s at 
the scale /ihard = Mq through the matching procedure. If Aqcd were of or- 
der of or even larger than Mqv (which is essentially the situation for cc) this 
would be all we could do using perturbation theory. For Aqcd smaller than 
Mqv, however, one can go further and also integrate out gluonic (and light 
quark) momenta of order AfQt>The resulting theory has been called "potential 
NRQCD" (PNRQCD) in Ref. and is characterized by the fact that its La- 
grangian contains spatially non-local four-fermion interactions which are noth- 
ing else than instantaneous (static) QQ potentials. The "short-distance" coef- 
ficients of the corresponding operators describe the gluonic (and light quark) 
effects from momenta of order Mqv and can be calculated perturbatively at 
the scale /i so ft = Mqv. To NNLO (i.e. including potentials suppressed by at 
most a 2 , oc s /Mq or 1/Mq relative to the Coulomb potential) the relevant QQ 
potentials read (a s = a s (/x so ft), Ca = 3, C F = 4/3, T = 1/2, ft = e 7 /z so f t , 
r = \r\) 
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where St and Sj are the top and antitop spin operators, L is the angular 
momentum operator and /3q,i are the one- and-two-loop beta-functions. The 
constants ai t 2 have been calculated in Refs. 0Q. V c is the Coulomb (static) 
potential and Vbf the Breit-Fermi potential known from positronium. V^ A 
is a purely non-Abelian potential generated through non-analytic terms in 
the one-loop vertex corrections to the Coulomb potential involving the triple 
gluon vertex. The remaining dynamical gluon (light quark) fields can only 
carry momenta of order Mqv 2 and describe radiation and retardation effects. 
If Aqcd < Mqv 2 one can show that those retardation effects are of NNNLO 
in a using arguments known from QED and taking into account how the 
gluon self coupling scales with v for gluonic momenta of order Mqv 2 . [This 
only works because the QQ pair is produced in a color singlet state!] This 
means that retardation effects (and the scale Mqv 2 ) can be ignored at NNLO 
and that the QQ dynamics can be described by a two-body positronium-likc 
Schrodinger equation equation of the form (E = y/a — 2Mq) 



Mo 



AM 3 Q 



V c (r) + V BF {f) + V NA {f) 



(5) 



containing the heavy quark kinetic energy up to NNLO and the instantaneous 
potentials (||)-(^|). In Eq. (^) Mq is defined as the pole mass. If Aqcd were of 
order or even larger than Mqv 2 , on the other hand, the coupling of the radi- 
ation gluon with the heavy quark would become of order one and retardation 
effects would be NNLO. In this case a perturbative calculation of <r thr would 
be impossible at NNLO because retardation effects would be non-perturbative. 
For this reason we have to make sure that the condition Mqv 2 > Aqcd is sat- 
isfied, if the NNLO expression for a thr derived in the following section shall be 
trusted. 



3 Calculation of the Total Cross Section 

For simplicity we only consider the photon mediated total cross section. The 
inclusion of the Z exchange is trivial for the vector current contributions. 
The contributions from the axial-vector current can be easily implemented at 
NNLO because the axial- vector current produces the QQ pair in a P-wave state 
which leads to a suppression oc v 2 relative to the vector current contribution. 
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We start from the fully covariant expression for the normalized total cross 
section (Rqq = <J thr /<r(e + e~ — > ^i + /i~)) 

#QQ(? 2 ) = ^M-i / ^e 4 ^(0|r^(x)/^(0)|0)] 

= iZ^i Ira [ _i ( |T #(g) ?"("?) | )] , (6) 

and expand the electromagnetic current which produces/annihilates the QQ 
pair with cm. energy yg^ in terms of 3 Si NRQCD currents up to dimension 
eight (i = 1,2,3) 

= b x (fta iX ) (g) - ^ (^V,(-§ B) 2 x) (?) + ..., 
Q 

H-q) = h (-«) - Jfj- (xV<(-i B) 2 ^) (-«) + . . . , (7) 

Q 

where the constants b\ and 62 are short-distance coefficients normalized to one 
at the Born level. Only the spatial components of the currents contribute 
at the NNLO level. Inserting expansion (0) back into Eq. (^J) leads to the 
nonrelativistic expansion of the NNLO cross section 

_ ^.2 

Al(E, jUsoft) Mfac) 

+ ..-, (8) 
where 

At = i(0\(ft* X ) (X j fi>)\0), (9) 
A 2 = \i{Q\{ft3x) (x^{-^D) 2 i;) +h.c.|0). (10) 

The cross section is expanded in terms of a sum of absorptive parts of non- 
relativistic current correlators, each of them multiplied by a short-distance 
coefficient. In fact, the right-hand side (RHS) of Eq_(||) just represents an ap- 
plication of the factorization formalism proposed inta. The second term on the 
RHS of Eq. m) is suppressed by v 2 , i.e. of NNLO. This can be seen explicitly 
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by using the equations of motion from the NRQCD Lagrangian, which relates 
the correlator A2 directly to Ai, 



A 2 = M Q EA 1 . (11) 

Relation (|ll|) has also been used to obtain the coefficient —4/3 in front of the 
second term on the RHS of Eq. (||). The nonrelativistic current correlators A\ t2 
contain the resummation of the singular terms mentioned previously. They 
depend on the renormalization scale /i S oft through the potentials (||)-(^). The 
constants C\ and C2 (which are also normalized to one at the Born level), on 
the other hand, describe short-distance effects and, therefore, depend on the 
hard scale /^hard- They only represent a simple power series in a s (where the 
coefficients contain numbers and logarithms of Mq, /if ac and /^hard) and do not 
contain any resummations in a s . At NNLO they have to be calculated up to 
order a 2 s because we count a s /v of order one for a perturbative nonrelativistic 
QQ system. 

The nonrelativistic correlators A1.2 are calculated by determining the 
Green function of the Schrodinger equation (||) where Vbf is evaluated for the 
3 S*i configuration. The NNLO relation between the correlator Ai and Green 
function reads 



Ai = 6N C 



lim G(r,f,E) 

\r I , \f' I — >0 



(12) 



Eq. ( |12[ ) can be quickly derived from the facts that G(r, r* ,E) describes the 
propagation of a quark-antiquark pair which is produced and annihilated at rel- 
ative distances \r\ and |r"|, respectively, and that the QQ pair is produced and 
annihilated through the electromagnetic current at zero distances. Therefore 
Ai must be proportional to lim^u^i^o G(f,f',E). The correct proportional- 
ity constant can then be determined by considering production of a free (i.e. 
a s = 0) QQ pair in the nonrelativistic limit. (In this case the Born cross sec- 
tion in full QCD can be easily compared to the imaginary part of the Green 
function of the free nonrelativistic Schrodinger equation.) The correlator A2 
is determined from Ai via relation ( |ll|) . We would like to emphasize that the 
zero-distance Green function on the RHS of Eqs. (|l2| ) contains UV divergences 
from the higher dimensional NNLO effects which have to be regularized. In 
the actual calculations carried out in Refs. Ln3 we have imposed the explicit 
short-distance cutoff /if ac , called factorization scale. This is the reason why 
the correlators also depend on fif ac . One way to solve Eq. (||) is to start 
from the well known Green function of the nonrelativistic Coulomb prob- 
lem and to incorporate all the higher order terms via first and second order 
Rayleigh-Schrodinger time-independent perturbation theory. 
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To determine the short-distance constant C\ up to 0(a 2 ) we can expand 
expression (g) in the (formal) limit a s <C v -C 1 (for ^t so f t = /Xhard) up to 
O(a^) and demand equality (i.e. match) to the total cross section obtained at 
the two-loop level in full QCD keeping terms up to NNLO in an expansion in 
v. In this limit fixed multi-loop perturbation theory (i.e. an expansion in a s ) 
as well as the nonrelativistic approximation (i.e. a subsequent expansion in v) 
are feasible. Because the full QCD cross section is independent of /Zf ac , C\ 
also depends on /Xf ac . We call this kind of matching calculation at the level of 
the final result "direct maching" . The consistency of the effective field theory 
approach ensures that C\ only contains contributions from momenta of order 
Mq and does not have any terms singular in v. Details of this calculation and 
references regarding the important calculations of the two-loop cross section 
in full QCD at NNLO in the velocity expansion can be found in Ref. B. A 
very economical method to calculate Feynman diagrams in full QCD in an 
expansion in v using dimensional regularization, the "Threshold Expansion" , 
has been developed in Ref.B. 

4 Applications 

4-1 Bottom Quark Mass from T Mesons 

Due to causality, derivatives of the electromagnetic current-current correlator 
with respect to q 2 at q 2 = are directly related to the total photon mediated 
cross section of bottom quark-antiquark production in e + e~ annihilation, 

= /;£*.<.>■ <13) 

Assuming global duality, P n can be either calculated from experimental data 
for the total cross section in e + e~ annihilation^ or theoretically using quan- 
tum chromodynamics (QCD). It is the basic idea of this sum rule to set the 
moments calculated from experimental data, P% x , equal to those determined 
theoretically from QCD, P^ h , and to use this relation to determine the bottom 
quark mass (and the strong coupling) by fitting theoretical and experimental 
moments for various values of n. At this point it is important to set the range of 
allowed values of n for which the moments P„ can be trusted using the NNLO 
cross section calculated in the preceding section. As mentioned, we have to 
make sure that Mb(v) 2 > Aqcd where (v) is the effective cm. bottom quark 
velocity in a particular moment. One can show that (v) ~ l/V"- f° r large n. 

b At the level of precision in this work the Z mediated cross section can be safely neglected. 
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Figure 1: (a) Allowed region in the e -a s plane for the unconstrained fit using the 

moments at NNLO. The gray shaded region represents the allowed region. The dots are 
points of minimal x 2 for a large number of models. Experimental errors arc included at the 
95% CL level and have been drawn as narrow gray ellipses centered at the dots of minimal 
X 2 . (b) Allowed M?° c values for a given value of a s . The gray shaded region corresponds 
to the allowed ranges for the NNLO analysis and the striped region for the NLO analysis. 
Experimental errors are included at the 95% CL level. It is illustrated how the allowed range 
for A/ 6 polc at NNLO is obtained if 0.114 < a s (M z ) < 0.122 is taken as an input. 



This means that n should be sufficiently smaller than 15 — 20. It is interest- 
ing that the same conclusion can be drawn from the Poggio-Quinn- Weinberg 
argument that the effective energy smearing range contained in the moments 
should be larger than Aqcd t3- On the other hand, n has to be large enough 
that the use of the cross section at threshold is justified because the energy 
regime close to threshold dominates P n only in this case. In our analysis □ we 
have taken the range 4 < n < 10. Larger values of n increase the danger of 
large systematic errors. To determine the allowed range for the bottom quark 
mass we have fitted P„ to P£ x for various sets of n's. It turns out that the 
theoretical errors are much larger than the experimental ones. The dominant 
theoretical errors come from the dependence of P„ on the scale /U so ft- We 
have combined experimental and theoretical errors by using the "scanning" 
method t3 in which a large number of statistical fits is carried out for various 
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Table 1: Recent determinations of bottom quark masses using QCD sum rules for the T 
mesons, m^mt) refers to the MS mass. NLO refers to analyses including corrections of 
order a a to the nonrelativistic limit and NNLO to analyses including corrections or order 
o? s , ot a v and v 2 . No order is indicated for Ref. 10 because the bound state poles have not 
been taken into account for P„ in that analysis. For entries where also values of a s are 
given a simultaneous fit for mass and QCD coupling has been carried out. In our analysis 
(Ref. 6) the errors are not written as Gaussian errors but as allowed ranges. 



authors 


order 


n 


Af h polc [GeV] 


m b (m b )[GeV] 


Rcf^ 


NLO 


4-10 


4.64 - 4.y2 
a s {M z ) = .086 - .132 






NNLO 


4-10 


4.74 - 4.87 
a a (M z ) = .096 - .124 


4.17 - 4.35 




NNLO 


4-10 


4.78 - 4.98 


4.16 - 4.33 


Ref^ 


NLO 


8-20 


4.827 ± .007 
a s (M z ) = .109 ± .001 




Ref.B 




8-20 


4.604 ± .014 

a s (M z ) = .U8±™1 


4.13 ±.06 
a s (M z ) = .120±;°J° 


Rcf.M 


NLO 


10 - 20 


4.75 ± .04 




Ref. El 


NNLO 


10 - 20 


4.78 ± .04 




Ref.y 


NNLO 


14 - 18 




4.20 ± .10 


Rcf.0 


NNLO 


8-12 


4.80 ± .06 





"reasonable" choices for the scales, each called "a model" . It is believed that 
this method represents a conservative way to combine experimental and large 
theoretical errors. Our result for a simultaneous fit for the bottom pole mass 
and a s is displayed in Fig. |l|a. In Fig. [j]b the result for the pole mass is dis- 
played if a s is taken as an input. It is conspicuous that the extracted values 
for the pole mass are quite different for both methods. This variation could be 
explained from the fact that the pole mass is nat.xLenned beyond an accuracy 
of Aqcd due to its strong infrared sensitivity. lIZHla. It therefore seems to be 
more advantageous to extract a short-distance mass like the MS mass. Taking 
into account the strong correlations between pole mass and a s and using the 
two-loop conversion formula between pole and MS mass we find very good 
agreement in the mass determination for both methods (see Tab. [I]). Using 
the moments at NLO we also found errors which..were much more conservative 
than the ones given in the NLO analysis of Ref. El. In Tab. [l] we give a compi- 
lation of all recent sum rule determinations of bottom quark masses based on 
experimental data from the T mesons. 
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4-2 Top Quark Pair Production Cross Section at Threshold 

The measurement of the tt production lineshape at threshold (y/s ~ 2M t ) is 
among the first tasks of the Next Linear Collider (NLC). Due to the large 
top quark width {T t (t -> Wb) w (G f /V2)M?/8tt w 1.5 GeV) which serves 
as a natural infrared cutoff and smearing mechanism one can in fact calcu- 
late the lineshape locally without imposing any additional smearing. To be 
more specific, the effective cm. velocity v Q g of the top quarks is of order 
y/E 2 + T 2 /M t which means that M t v 2 s is larger than Aqcd for any nonrela- 
tivistic energy. Physically this cutoff mechanism arises because the top quarks 
decay weakly before hadronization effects can set in. This also leads to the 
phenomenon that there are no individual narrow toponium resonances because 
the latter are so broad that they are almost completely smeared outlia In other 
words, there will be no toponium spectroscopy as we know it from charmo- 
nia or bottomonia. However, the large decay width allows for remarkably 
accurate measurements of top quark properties without the diluting effects of 
hadronization. It is believed that measurements at the tt threshold at the 
NLC will provide the most precise determination of the top quark mass. In 
view of the potentially sizeable NNLO relativistic corrections it is mandatory 
to calculate and control the NNLO effects. Our calculation of the NNLO cross 
section has been designed for stable quarks and can, strictly speaking, not be 
used as the NNLO cross section for tt^\ Nevertheless we can illustrate the 
impact of the NNLO corrections by using the naive replacement E — ► E + iTt 
in Eq. (^) where T t is the free top quark width—. This prescription has been 
proven to be correct in the nonrelativistic limit and should be sufficient to 
give us the bulk of the NNLO relativistic corrections to the ti cross section 
at threshold. In Fig. | the LO (dotted lines), NLO (dashed lines) and NNLO 
(solid lines) normalized cross sections are plotted versus E = \fq^ — 2Mf olc in 
the range -5 GeV< E < 5 GeV for A/ t pole = 175 GeV, a s {M z ) = 0.118 and 
F t = 1.43 GeV. For the scales the choices /i so f t = 50 (upper lines), 75 (middle 
lines) and 100 GeV (lower lines) and /ihard = Mfac = M t have been made and 
two-loop running of the strong coupling has been used. It is evident that the 
NNLO corrections are large. The behavior of the NNLO corrections clearly 
indicates that the convergence of the perturbative series for the tt cross section 
is much.. worse than expected from the general arguments given by Fadin and 
Khozelij. Further examinations of the NNLO cross section including a study of 
the impact of the use of a mass other than the pole one and a proper treatment 
of the top width are mandatory. 



c As already mentioned, we also have not yet included the Z exchange contributions. 
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1.2 




E(GeV) 



Figure 2: The total normalized photon-mediated ti cross section R t j at LO (dotted lines), 
NLO (dashed lines) and NNLO (solid lines) in the nonrelativistic expansion for the soft 
scales /i so f t = 50 (upper), 75 (middle) and 100 GeV (lower). The plot has been calculated 
in Ref. 7. Similar results have been obtained in Rcf. 20 and 21. 

5 Conclusions 

During the last year there has been significant progress in the understanding 
of perturbative heavy quark-antiquark systems in the kinematic regime close 
to threshold. Using the concepts of effective field theories our knowledge has 
increased at the conceptual level and a number of previously unknown NNLO 
corrections have been determined. In this talk I have reviewed the ideas in- 
volved to perturbatively calculate and then apply the NNLO corrections of the 
total QQ production cross section at threshold. 
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